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ABSTRACT 


Advances in the field of integer programming have 
recently been made by Young and Glover, who have developed 
primal all integer integer programming algorithms. It is : 
the purpose of this paper to gain computational experience 
and suggest techniques to improve the efficiency of the 


algorithm. A first phase procedure is developed to obtain 


an initial feasible integer solution. 
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I Introduction 
Advances in the field of integer programming have 
[10] [5] 


recently been made by Young and Glover » who have 


developed primal all inter-integer programming algorithms. 


A primal algorithm is one that proceeds through a 
series of primal feasible bases towards an optimal solution 
of the problem. The two principle advantages of a primal 
method are: 

1) at any point during the solution, the problem may 
be stopped and a feasible solution exists that may be used 
as a "good" approximation: 

2) a feasible solution arrived at by any method may 
be used as a starting point for a primal algorithm. The 
most famous example of a primal programming algorithm is 


Dantzig's Simplex Algorithm. '?-34 


[9] 


Until Young developed his first primal algorithm 
there were no finite primal algorithms available in the 
literature; however, there was one privately attributed 
to Gomory. ‘Algorithms developed prior to this time were 
either cutting plane or enumeration types. Young's first 
algorithm was extremely complex and difficult to implement 
for any practical purpose. Subsequently Young simplified 


his algorithm!101, 


co 
The objective here is to gain computational experi- 


ence with the new algorithms of Young and Glover.and to 


suggest techniques that will improve their efficiency. 





The two algorithms are not significantly different, thus 
they will be considered as one where the minor differences 
will be noted. 

The algorithm contains nondeterministic rules which . 
will be investigated to determine their effect on the pro- 
gress of the algorithm. The rule for the selection of a 
pivot column in the Simplex Algorithm for solving linear 
programs is an example of a nondeterministic rule. In 
its most general form the rule is to select any column that 
has a negative coefficient in the objective function. The 
user has to determine which column he will select as a 
pivot column from a set of columns that have negative co- 
efficients in the objective function. 

The algorithms require a knowledge of a feasible in- 
teger solution. In a problem this may be difficult to Si 


obtain . A first phase procedure is developed to obtain an 


initial feasible integer solution, 








II Description of the Algorithm 
The basic problem under consideration is: 
find Xa j=1,2,...,n that maximizes Xo 


where 


Tal 
n 
Xs + ) a5; Ss os i=192,<.%,m 
gall 
(1) 
a, = Se Se weigh 


ss 20 and integer JEly2e.s eet 


Without loss of generality it can be assumed that all 
ais are integers and all as 2 0.° If not, appropriate 
techniques are available to insure this. 

In this section the procedures to solve (1) are de- 
scribed and set forth. The method is contained in a primal 


algorithm that relies heavily on Dantzig's Simplex Algo- 


[2,3] and Gomory's All Integer Integer Programming 
(7] 


rithm 
Algorithm 
Essentially the algorithm selects a pivot column and 


Ce] 


a source row which is used to generate a Gomory cut The 


cut equation is then adjoined to the problem and used as 


the pivot row for a change of basis. Each pivot does not 





necessarily result in an increase in the objective function 
and more restrictive rules are necessary for the selection 
of a pivot column and source row to insure that the proce- 


dure is finite. 


Prior to proceeding to the description of the algo- 
rithm it is necessary to define the following terms and no- 


tational conventions. 


Let 
v = index of pivot column 
Ss = index of source row 
J = set of indicies of non basic variables 
8, = . : ae ee ay | oe 
j min (A g/g) for a So  £6F° 9 and. 


[t] = greatest integer less than or equal to t. 
S(3) = set of row indices i such that 
o<la, /a,i| = 8, for 7 ih Js 
io ij 7 
A. = jth column of coefficients for j in J. 


The cut equation is defined as 


de, / aay | , Ti tes Laso 4 acy (2) 
J¢€T 


where u = slack variable associated with the cut 
equation 
A transition cycle is a pivot that results in an in- 


crease in the objective function. A stationary cycle is 


a pivot that does not increase the objective function. 





The reference constraint is a constraint of the form 
. x. *< 
y ay 5X5 aro (3) 
J¢T 
that is adjoined to the problem to guide the progress of 
the algorithm through a sequence of stationary cycles. 
This is necessary to insure that the number of interven- 
ing stationary cycles between two transition cycles is 
finite which is necessary to insure that the algorithm will 
teminate in a finite number of steps. 


The reference vector is defined as 


as 
ry 
for all j such that a5 > 0. A, is not defined for 


a_.' S 0. 
ry 


The exact form of the reference constraint is not spe- 
cified and is one of the choices of the user. The refer- 


ence constraint must satisfy the following. 


1) A a implies a > o~ 
j ” rj 
2) <0 impli A. >A (5) 
ek es : 
rj impli = 
a * The symbol < means lexicographically less than. A vector 


x is said to be lexicographically less than a vector y 


if the first non zero component of the vector (x-y) is 


negative. 





Young suggests using (3) as a reference constraint 
where all aps = 1 and Be is a positive constant selected 
to be large enough so that no feasible solutions are ex- 
cluded from the problem. Glover suggests using a surro- = 
gate constraint as the reference constraint. 

A surrogate constraint is “an inequality implied by 
the constraints of an integer program, and designed to 
capture useful information that cannot be extracted from 
the parent constraints individually but is nevertheless a 


[5] 


consequence of their conjunction". In general a surro- 
gate constraint is of the form of equation 8 below and is 
obtained from a linear combination of the constraints of 


the original problem. 


If the original problem (1) is written in the form: 


; Ps T 
find x = (X)+Xoe ome x.) that 
maximizes Xo = cx (6) 
where Ax <b 
and x 20 and integer 
when Cr (854 7899° eee a Aon)» A — aly a = diyj2oe ae M, 


} = 1ly2,ce0 FB 


a a 4 
LO’ 3207 i  hmOr ae 


dropped, the dual of this problem is then: ‘ 


and b= (a If the integer constraint is 
find w= (W)+Wo, aks wi) 


that minimizes w, = Wy i 


0 
(7) 


where wA 2c 
w 20 





[4] 


The reference constraint as suggested by Glover is 


then defined as 


n 
a3%5 <a. (8) 
. j=1 
where a. = Ww*aA 
rj J 
a = w*b 
ro 
w* = a feasible solution to the dual. 


An appropriate source row selection rule is defined 
as any rule which guarantees, for any row k, that atfinite 
intervals a tableau occurs in which a,,, $ a¥- The rule 
is necessary to insure that a transition cycle will occur 
at finite intervals. 
The following rules completely determine the algorithm:. 
1) Adjoin a reference constraint to the problem. 
2) Check for optimality. The current tableau is opti- 
mal if 554 20 for jedg or if Boy ary > -1/a..- 
3) Determine the pivot column v, such that AY <3, 
for j # v and a5” 0. 
4) Select a source row, s, from S(j) by an appropri- 
ate source row selection procedure. 
5) Adjoin to.the tableau the Gomory cut (2) derived 
from the source row. 
6) Carry out the usual simplex change of basis with 


vas the pivot column and the Gomory cut as the 


pivot row. 





7) Delete the pivot column and pivot row from the 
system. 


8) Go to 2. 


Steps 5,6, and 7 may be combined to one step as: 
5a) Execute the change of basis as follows: 


1) a.. = /- aa for j =v 


NO 
— 
Il 


ij ai5 7 a,5/ey] e&,. tery #v (9) 


Step 5a is equivalent to steps 5,6, and 7. 


The above procedure always selects the column with 
the lexicographic minimum A, as the pivot column. It may 
be desirable to execute a transition cycle when possible 
in order to immediately increase the objective value and 
possibly decrease the number of iterations. In this case, 


Rules 3 and 8 above are replaced by: 


3a. If for some j in J and A654 <0 5 21, designate 
that column as the pivot column. If none exists 


use rule 3 above. 


8a. Go to l if a transition cycle has been executed. 


Go to 2 otherwise. 


This procedure requires that the reference constraint 
be revised after each transition cycle since the problem is 
essentially restarted as if it were a new problem. Because 


of this, computationally it would be advantageous to have an 


easily determinable reference constraint. 





Optimality is achieved if dual feasibility occurs, 

is i 2 i mm = ‘ 
that is if all ae 0, or if ae 1/a., The second 
criteria occurs if no positive integer values exist for the 
non basic variables that will satisfy the reference con- 
straint and improve the solution. A complete proof is con- 


[4]. 


tained in 


Rules 3 and 4 are the essential rules necessary to in- 


sure that the algorithm is finite. Basically, these two 


(x) the reference vector for the ~ 


iteration, is lexicographically greater than Hy rls 


rules insure that A, 
the 
reference vector for the previous iteration. Since the 
first component of AL, is a and AL, A. for j # v and 


. j in J, for each iteration, then when A> 0, the problem 


is at optimality. 


15 





III Discussion of Problem 

Within the algorithm there are three non deterministic 
rules, the form of the reference constraint, the source row 
selection rule, and whether to execute a transition cycle 
as soon as possible or to select the column with the lexi- 
cographically minimum reference vector. 

The reference constraint must satisfy the criteria of 
(5.1) and (5.2). One such constraint is that of (3). 
Denote this by Reference Constraint One (RC1). The cri- 
tical point is determining aoo° Any large positive inte- 
ger can be used that does not exclude any feasible integer 
solution. .It is possible to determine the upper bound on 
the summation of the variables by solving the linear program 
of (6) where the integer constraint is dropped, c=(1,1,...1) 
and all other variables are as defined in (6). 

A second form of the reference constraint is that of 
(8), a surrogate constraint, as suggested by Glover!4), This 
will be denoted by Reference Constraint Two (RC2). 

The final form of the reference constraint considered 
is: 

cx s Fo (10) 

where cx is the original objective function and my is some 
large integer constraint such that no feasible integer 
solutions are excluded from the problem. Denote this form 
as Reference Constraint Three (RC3). As in RCl, c, can be - 


obtained from the solution of the linear program (6) with- 


out the integer constraints. 


16 





If a transition cycle is executed out of order, that 
is when first possible rather than always selecting that 
column with the lexicographically minimum reference vector, 
then the reference constraint must be revised after each 
transition cycle. Because of this RC2 was not considered 
when transition cycles were executed out of order. It did 
not appear to be a reasonable procedure to solve the dual 
problem after each transition cycle. 

To meet the criteria of the definition of an appro- 
priate source row selection rule, it is necessary to insure 
that for any row i only a finite number of cycles will 
exist such that a._ 2 a... In order for this to occur, 


aN LO 


once a row has been selected such that asy 2 Aso? it is 
necessary to repeatedly select this row. until it no longer 
qualifies as a source ro\W. 

The source row selection rules (SRR) that were used 
are 


SRR1) a) Select that row in S(v) with the smallest 


row index number. 


b) Repeatedly select this row until it no 


longer qualifies. 


c) go toa 
SRR2) a) Generate a sequential list of row indicies. 
b) Select the first index of the list that is 

in S(v). 


c) same as SRRIb. 
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d) Continue down the list until a row is 


found that is in S(v), go toc. 
e) If none found before end of list go to b. 


SRR3) a) Generate the vector G, for each row in 


S(v) where 
g(j) = ai5/Aivl I = ligetuw B (11) 


b) Select source row s such that G. 5 G; 


for all iin S(v) i¢s 
c) same as SRR1b 
dad) Go toa 


SRR1 is used because of the computational simplicity. 


[10] SRR3 is 


SRR2 is one of the rules Young suggested 
used to increase the values of the coefficients of the ob- oe 
jective function. As a justification for this, consider 


rule 5a of the algorithm, which defines the pivoting process 


for the objective function row, 


«| Si i #Vv (12) 


‘ a fa | a 
oF oF sj’ svi “ov 


It can be seen that B54 will be changed by an amount 
he ; ‘ # Ae te 
$3) aoy in the next tableau. Since aov 0, Te we 
desired to have G. as large as possible (in some sense), in 


order that Boj will increase towards the optimality criteria 


of alla... 20. . 
oO} 





With the above choices, fifteen different combinations 


were considered on twelve test pcubiae eo The next sec- 


tion describes the problems and the results. 
The algorithm is demonstrated by means of a sample 


problem in Appendix 2. 
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IV Computational Results 

The results of the various choices discussed in the 
last section will be presented here. The characteristics 
of the twelve test problems are summarized in Figure l. 
Problems 7 and 8 are identical to 5 and 6 respectively with 
two constraints omitted. The constraints are 0-1 bounds 
on two variables. For problems 9, 10, and 11, RC2, the re- 
ference constraint derived from a surrogate constraint, was 
exactly the same as RC3, the reference constraint that uses 
the objective function as a constraint. 

Figure 2 is the summary of the results of the problems. 
RC3 with transition cycles executed as soon as possible was 
considered and a few problems were run but the results were 
so poor that it was not pursued any further. The reason 
for the poor results was that the optimality criteria of » 
ee fe = -1/a.. was never met. This resulted because 
Bo Fry was always equal to -1 and ~1/a,. was a negative 
fraction. This made it necessary to achieve dual feasibili- 
ty, before the problem could be terminated, which consider- 
ably increased the number of interations. 

The computational procedures are straight forward and 
are completely described by the rules in Section three. 
One variation considered was that of obtaining an initial 
basic feasible solution to the problem when one was not 


readily available. : 
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Problem Number Number Type* 
of of of 
Number Constraints Variables Constraint Denominator ** 


« 1 4 5 all LE 183 
2 4 5 all LE 258 
. 3 4 5 all LE 320 
4 4 5 all LE 205 
5 6 5 all LE 19400 
6 6 5 all LE 32000 
7 4 5 all LE 19400 
8 4 5 all LE 32000 
9 6 6 all LE 2000 
10 "] ¥ all GE 32 
LL 7 7 all GE 32 
12 2 4 all GE 72 
Characteristics of the Test Problems 
: * LE indicates that the constraints are less than or 


equal to inequalities. 
GT indicates that the constraints are greater than 


or equal to inequalities. 


** The denominator is the determinate of the coeffi- 
cients of the basic variables when the problem is 
solved as a linear program. It gives a general indi- 
cation as to the difficulty of finding an integer so- 


lution. 


Figure One 
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ce cee ere re ces me se ne ee ee es is ee ere ee we ee ae we ee ee a i ee es ee ee i 


** 


67 77 854 401 387 168 316 287 391 303 


59 24 422 1021 355 ALS 569 188 B15 419 
537 171 479 1124 480 339 1202 520 29 29 
5 5 399 24.12 325 623 691 48 39 46 
90k 26533 90k 90k, 31514 90k 90k 90k 46876 47669 
90k 90822 90k 90k 90k 90k 90k 90k 234 234 
90k 90k 90k 90k 203 90k 90k 90k 54075 53867 
90k 90k 90k 90k 90k 90k 90k 90k 234 234 
600 240 12289 19757 2548 12289 18757 2548 3684 4936 
18 18 31 20 24 39 20 24 18 33 
43 25 34 120 39 34 120 39 19 22 
11798 3822 226 Lé65 116 39 39 7) 21 451 


Number of “Iterations to Solve the Test Problems. 


Pivot column always the column with the lexicographically minimum 


reference vector. 
Transition cycles executed when possible. 


A k indicates that the problem had not terminated in the indicated 
number of iterations (x10). 


Figure Two 
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This occurred in problems 10, 10, and 12. A Phase I tech- 
nique analogous to that of the Simplex Agorithm was used. 
Rule 5a was used in place of rules 5, 6, and 7 of the algo- 
rithm for the computations. This made it necessary to 
employ a special procedure to eliminate artificial variables 
from the problem during Phase I. At various points during 


Phase I, equation exist of the form 
“u; + y. = 0 (13) 


where uy is a slack variable generated from a cut and Yi is 
the co artificial variable. When equations of this form 
appear during Phase I they are dropped and the a column 
is deleted from the problem. This is done because the 
artificial variable, Yj equals zero and the equation implies 
u, equals zero. It is then eliminated from all other equa- 
tions in the problem. One further special procedure was 
employed during Phase I. When w is reduced to zero, it is 
possible that some of the artificial variables are still in 
the problem. They must be deleted before proceeding to 
Phase II. To accomplish this a pivot is carried out on any 
element that will give the form of (13). This can best be 
accomplished by pivoting on an element where some V3 is 

the basic variable and the coefficient is equal to one, 

Any column may be selected as the pivot column if the pivot 
is carried out in an equation that has an articial as a 
basic variable. This is possible because a... = © for all 


rows i that have artifical variables as basic variables 
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(when w= 0). This special procedure is carried out until 


equationsof the form of 13 are obtained and are then de- 


leted from the system. When all artificals have been re- 


moved an initial basic feasible solution has been found 


and the problem can be solved by the basic algorithm. 





V Discussion of Results 

From the summary of results in Figure 2 a number of 
conclusions can be made. It appears that the best proce- 
dure investigated is to execute a transition cycle when 
first possible in order to reduce the number of iterations. 
The best reference constraint, of those investigated, was 
RCl as originally suggested by young/202 | This constraint 
produced the least number of iterations in general. It 
also was the easiest to obtain and use in terms of compu- 
tational simplicity. This result was somewhat unexpected. 


[4] 


It was initially believed that RC2, the surrogate con- 
straint, would yield better results due to the nature of a 
surrogate constraint. This advantage however was complete- 
ly overshadowed by the advantages gained by executing a 
transition cycle as soon as possible. RC2 was not investi- 
gated when a transition cycle is executed out of order. It 
does not appear to be a reasonable procedure to solve the 
dual problem after each transition cycle, as would be re- 
quired, 

The results of the various choices of source row 
selection rules were not as clear as that of the reference 
constraints. The only rule that significantly decreased 
the number of iterations was SRR3 when used in conjunction 
with RCl, and when transition cycles are executed out of 


order. For the other cases it appears that SRR-3 is 


slightly better than the other two rules. 
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All of the problems tested here were of moderate size 
yet the number of iterations necessary to reach optimality 
was usually quite large. On larger problems the number of 
iterations may be so large that the algorithm would be of 
no practical use, even when the most ingenious and effi- 
cient techniques are used. 

One of the most promising aspects of the algorithm is 
that it may open the way for a variety of composite or 
specialized algorithms to fit the needs of specialized users. 

Among the areas that need future investigation, the 
following appear to be most significant: 

1) More work is necessary to reduce the number of 
iterations necessary to reach optimality. On a number of 
the test problems, optimality was achieved in relatively few 
iterations, but the criteria was not met for many more ite- ‘ 
rations. On test problems 1 and 2 optimality was achieved 
after one iteration but even in the best cases at least 20 
or more iterations were necessary to recognize optimality. 
Glover has made an outstanding contribution toward this end 


[4] 


with his modified optimality criteria as used here. 


2) Another area, which requires future work is devis- 
ing a method that requires less computer storage and more 
importantly, fewer calculations. The basic tableau re- 
quires a space of (m+ n+ 2) x (n + 1). It would not take 
a very large problem to exhaust the memory of a present day > 


computer. If the necessary number of calculations could be 
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decreased, a larger number of iterations would be tolerable. 
In test problems 5 and 6 it took approximately 10 minutes 
of execution time to carry out 90,000 iterations. 

3) Further work is also necessary to determine an 
initial basic feasible solution. It has been shown here 
that a Phase I technique can be used successfully. 

The algorithm is a major step forward in the field of 
integer programming. In theory, it provides a finite algo- 
rithm that will solve any integer programming problem. In 
practice the number of iteration necessary to reach opti- 
mality can be large, especially on the more difficult pro- 
blems if the choices for the non-deterministic rules are 


not made carefully. 
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Appendix A 

The following problem will be worked as an example of 
the algorithm. 

Find x =(X)7XorX3eXqeXe) that 


. aximizes xX_=x,+x,+x 
maximi ok 3 tkytks 


where 2x, +3X5 +x,4+2x,+2x, +) = 18 

3X, +2x5+2x,+ X,+2x, +Y5 = 15 

~6x, SE Xs +Y3 = 0 

-7X5 +X) +Y, = 0 

~Xy xy = 0 

“Xo +X5 = 0 

~X3 +X, = 0 

“Xy +X, = 0 

. “Xp +X, = 0 
The reference constraint is 
X +X, tk, 4K, tk, +t. = 11 

where a = 11 was determined by solving the linear program 


ro 


as described in Section 3. The first tableau is: 


X, Xo X3 X4 Xp, b 
2 3 1 2 2 18 
3 2 2 1 2 s 
-6 0 a 0 0 ws 
QO -7 ) 1 0 @) 
-l -7 0 ) 0 ) 
. 6 =-—k po*' OB 0 
0 Oo -l1 0 0 0 
0 6) 0 -l 0 0 
: : : Fr ot reference constraint 
) Oo -1 -1 -1 0 objective function 


where the identity matrix has been omitted. 
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A transition cycle can be executed immediately by se- 
lecting column 5 as the pivot column. The minimum ratio 
0. = 7% and the source row is row 2. A pivot is carried 


out using rule 5a of the algorithm. The next tableau is: 


0 1 -1 2 <2 4 

2 0 0 1. ee si 
=5 0 1 0 0 0 

0 -7 0 1 0 0 

=1 0 0 0 0 0 

De at 0 0 0 0 Tableau 1 
0 0 <1 0 0 0 

0 0 0 = =1 0 0 

1 i i 0 1 7 

0 0 0 2 \) = 4 RC 
ah 1 O Sal nh g OF 


Since a transition cycle was executed, the reference 
constraint must be modified to: 


a. 1 d i 1 jal 





The problem is acutally at optimality now but it does 
not meet the criteria for optimality yet. The pivot column 
for the next iteration is column 4 (the only column eligibld. 
6, = 0 and the only row in $(4) is row 4 A stationary 


cycle must be executed. The next tableau is: 


0 15 -1l -2 -2 4 

1 7 Oo -1 -2 ai 
—6 (0) 1 0 0 ) 
-l 0 0 0 0 0 

Oo <1 0 @) 0 0) Tableau 2 

0 Oo -l 0 (@) 0 

QO -7 0 Bl 0 0 

By 1 IF 0 x 7 

1 6 ol ee 1 me 
1 -6 0 x ul 7 

Again only one column is available as the pivot column, . 

column 2. 6, = 1/7 and S(2) = (1,2). Using the criteria of 


2 
SRR3 the two G vectors are: 





G(1) (0,1,-1,-1-1,0) 


G(2) 


Il 


(0,1,0,-1,-1,0) 


therefore row 2 is selected as the source row. Again a 


stationary cycle must be executed. The next tableau is: 


1 
a 
a 

! 


13 13 
6 5 
i 0 
0 0 
=] =]. 
0 0 
7 
2 
9 
5 


FPREROOOFOFO 


Tableau 3 
=6 s 
1 
7 
=5 = 


e 
NPNOOO0OOFr A 


NOrAIOrOO 
l 
OrFRrFOROCOOCOF 


Two columns are eligible as pivot columns, 4 and 5. 


The first components of the reference vectors are: 


A(4) 


(-5/7,...-) 
- A(5) = (=5/9,.e.-) 


Since A(4) x A(5), column 4 is selected as the pivot column. 
6, = 1/6 and S(4) = (1,2). Row 2 was the source row for 
the last tableau so it will be selected as source row as 


long as it is in S(v). The next tableau is: 


0 aly) -l -13 13 4 
alt 5 0 -6 5 1 
-6 0 1 0 0 0 
0 -2 0 1 (6) 0 
-1 0 0 0 (0) 0 
8) -1 0 at ~1 6) Tableau 4 
0 0 -1 0 fe) 0 
0) -5 ) 6 ~7 6) 
al al 1 ~l 2 2 
1 6 1 -7 9 AEA 
al -4 0 5 -5 vj 


All aspects of the algorithm have been covered in the 


first four tableaus. The following tableaus will be pre- 
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The pivot element will be 


sented without further comment, 
listed below each tableau. 


The pivot element for the fourth 


tableau is (2,2) 


Tableau 6 


Tableau 5 





tat 2 BOA AQF ci & 


PES SONS er wh 


oy one Ha tae ™ 


ip ele @ Bird Satire 


Gets ea Peer 


Toi 20.990 OPE 


Honoo0ononhto 


(2,2) 


Pivot Element = 


(2,4) 


Pivot Element 


Tableau 8 


Tableau 7 


Tif See oo har 


NOONDOONHMH 


=i deel datiaa cethie l 


a ate Gee 2A Geis Raine 


SS al EON AED eS) ira Rh 


(4,4) 


Pivot Element 


Pivot Element = (4,5) 


Tableau 10 


Tableau 9 


SAGE S © 4 kent 
OoTtTooOonoyrTnrdhr: 
N 1 ain 
! I 
On Ee | ee 


cniasiliee tharuity’ alah Ai 
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iy eee ae ae 


4 
0 
0 
0 

-2 
(8) 


-11 20 
-2 
0 
-1 
0 
i 
0 
8 
-4 





(1,4) 


Pivot Element 


(10,5) 


Pivot Element 
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Tableau 12 


Tableau 11 


PaO eS. 22> a) 
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1 ow 7 ri! 
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Bad NOD eee Se ealied rd 


(1,4) 


Pivot Element 


(1,5) 


Pivot Element 


Tableau 14 


Tableau 13 
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el) ed SS ett 
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Pivot Element = (4,5) 


(4,3) 


Pivot Element = 


Tableau 16 


Tableau 15 
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-- 
Ww 
. 
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eo 
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Pivot Element 


(10,3) 


Pivot Element 
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Tableau 18 


Tableau 17 


Sep 2e2n tii 
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= 


Sa ad et 


Pivot Element = (1,3) 


(2,5) 


Pivot Element 


Tableau 20 


Tableau 19 
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| te 
1 
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til il 
Do | 
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tag ey 
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Pivot Element 


(10,2) 


Pivot Element 


(4,4) 


Tableau 22 


Tableau 21 


Thee Oe a2 err 
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St OO. DO fsal Pe 
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(3,2) 


Pivot Element 


Pivot Element = (3,4) 
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Tableau 24 


Tableau 23 
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OHFAGCFOOOFnANO 


i dl a 
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(10,1) 


Pivot Element = 


(1,4) 


Pivot Element 


Tableau 26 


Tableau 25 


BE PIA ED EL rel Fe 


Se OF OD rats © 


MWHAOUWUHANONKR AH 
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Pivot Element = (6,2) 


Pivot Element = (6,4) 


Tableau 28 


Tableau 27 


PASS Oi SS mel 
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OHAOFA CONF FARANO 


Final Tableau 


(4,4) 


Pivot Element = 
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Tableau 28 meets the optimality criteria of 
a- i - Pr 
B sey ay 1/ao? that is -2/23>-1/11, although the problem 
is not dual feasible. 


The optimal values of the variables are: 


xy = 7 
x, = 0 
X5 = 0 
x3 = 0 
Xx, = 0 
X_ = ) 


and the values of the slack variables are: 





: 
z 
— 
te 
= 
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= 
Wn 
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